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Abstract — The probability hypothesis density (PHD) and multi- 
target multi-Bernoulli (MeMBer) filters are two leading algo- 
rithms that have emerged from random finite sets (RFS). In this 
paper we study a method which combines these two approaches. 
Our work Is motivated by a sister paper, which proves that the 
full Bayes RFS filter naturally Incorporates a Poisson component 
representing targets that have never been detected, and a linear 
combination of multi-Bernoulli components representing targets 
under track. Here we demonstrate the benefit (In speed of 
track Initiation) that maintenance of a Poisson component of 
undetected targets provides. Subsequently, we propose a method 
of recycling, which projects Bernoulli components with a low 
probability of existence onto the Poisson component (as opposed 
to deleting them). We show that this allows us to achieve similar 
tracking performance using a fraction of the number of Bernoulli 
components {i.e., tracks). 

I. Introduction 

In recent years, random finite sets (RFS) fl] has emerged 
as coherent approach for inference in problems commonly 
encountered in tracking, involving unlabelled measurements 
of a set of objects with unknown cardinality. The probability 
hypothesis density (PHD) |2| and multi-target multi-Bernoulli 
(MeMBer) |1 3| filters have been shown to be effective in 
a variety of tracking problems. This paper demonstrates the 
benefits associated with combining these two structures in a 
hybrid Poisson-MeMBer tracker. The work is motivated by the 
derivation in the sister paper |4|5| , which shows that, under the 
common assumption that the target birth process is Poisson, 
the full Bayes RFS filter consists of both a Poisson component, 
and a linear combination of multi-Bernoulli distributions, 
where the Poisson component represents the distribution of 
targets that have never been detected. We show the benefit 
of incorporating both of these components in their natural 
roles, and subsequently show the benefit of extending the 
Poisson component to represent tracks {i.e., multi-Bernoulli 
components) with a low probability of existence. 

At first glance, it may seem unusual to be maintaining 
a distribution of targets that have never been detected. The 
reasons for doing so may be best understood in the context 
of agile sensors such as phased array radars. Intuitively, if 
the radar has not observed a region of space for a long 
period, it is more likely that there will be new targets in 
that portion of space awaiting first detection. Conversely, if 
a region has been observed very recently, we expect that it is 
comparatively unlikely that undetected targets will be present. 
The mathematical model behind this intuition is that targets 
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arrive (and depart) at a constant rate regardless of whether 
they are detected. We thus need to capture this arrival and 
departure of targets even in portions of space we are not 
presently observing in order to be able to predict the expected 
number of new targets in the area when it is next observed. 

Practically, the distribution of undetected targets is useful 
for two purposes. Firstly, it can be used for sensor scheduling, 
steering an agile sensor in order to minimise the expected 
number of undetected targets. This has been studied recently 
in fS], and is not considered in this paper Secondly, it has been 
shown |4| that the Bayes RFS filter naturally incorporates the 
distribution of undetected targets into its calculation of the 
probability that any detection represents a new target. Thus, 
maintaining a distribution of undetected targets allows the 
Bayes RFS filter to achieve the best possible track initiation 
performance, since it utilises the information to anticipate 
the relative likelihood of a measurement representing a false 
alarm, a target currently under track, or a new target. 

The multi-Bernoulli representation has significant advan- 
tages compared to the PHD as it is able to accumulate high 
confidence in the existence of targets. The major disadvantage 
is the requirement to maintain a large number of tracks 
with a small probability of existence in order to achieve 
adequate track initiation performance. The hybrid Poisson- 
multi-Bernoulli representation proposed in this paper moti- 
vates a true hybrid PHD-MeMBer algorithm which uses the 
Poisson component not only for undetected targets, but also for 
tracks with low probability of existence. This hybrid approach 
is shown to result in improved track initiation performance 
over the conventional method while maintaining only a small 
fraction of the number of tracks. A similar concept was 
proposed recently in 171; we discuss the relationship with this 
work in Section |Vl 

A. Outline and contributions 

Section III] describes the necessary results from the sister 
paper ||4|(which shows that the Bayes RFS filter naturally 
incorporates both Poisson and multi-Bernoulli components), 
concentrating on the role of the Poisson component and its 
interaction with the multi-BernoulU component. The remaining 
sections provide new contributions including: 

• Demonstration of the benefit of the hybrid Poisson- 
multi-BernouUi representation, in both uniform, station- 



ary problems (Section III-A i, and in problems involving 
agile sensors (Section III-B i 



Principled development of a hybrid PHD-MeMBer al- 
gorithm, providing a method for projecting a subset 
of multi-Bernoulli tracks onto the Poisson component, 
and an analysis of the distortion that this causes (in 
turn, providing insight into the recommended existence 



threshold to use in the projection) (Section IV-Ai 



Demonstration of the practical benefit of the hybrid PHD- 



MeMBer approach (Section IV-B i 



The relationship of the proposed method to prior work is 
described in Section [V] 

II. Background 

A. Assumptions and notation 

The assumptions and notation include: 

• The state of objects is denoted by x G A" (e.g., position 
and velocity in two dimensions) 

• Targets arrive according to a non-homogeneous Poisson 
point process (PPP) with intensity X^{x), independent of 
existing targets 

• Targets depart according to independent, identically dis- 
tributed (iid) Markovian processes; the survival probabil- 
ity in state x is P^{x) 

• Motion for each target is governed by a Markovian 
process, independent of all other targets; the single- 
target transition probability density function (PDF) is 
ft\t-i{x\x') 

• Each target may give rise to at most one measurement; 
the probability of detection in state x at time t is P^{x) 

• Each measurement z ^ Z {e.g., position in two dimen- 
sions, or range and azimuth) is the result of at most one 
target 

• False alarms arrive at time t according to a non- 
homogeneous PPP with intensity X^^{z), independent of 
targets and target-related measurements 

• Each target-derived measurement is independent of all 
other targets and measurements conditioned its corre- 
sponding target; the single target measurement Ukelihood 
is ft{z\x) 

We denote by Zt = {z}, . . . ,zl^'} the measurement set at 
time i, and by Z* = {Zi, . . . , Zt] the measurement history 
up to and including time t. The multiple sensor case may 
be addressed by performing update steps for each sensor 
sequentially in between prediction steps. 

B. Review of results from sister paper 

The sister paper |4| provides a form of conjugate prior for 
the tracking problem under the above assumptions. The form 
involves the union of an independent Poisson process and a 
linear combination of multi-BernoulU distributions. A practical 
algorithm is obtained by approximating the linear combination 
of multi-Bernoulli distributions as being multi -Bernoulli. The 
components of this form are: 

• a PPP representing undetected targets, with intensity 
\\t,{x), and 



a series of BernoulU tracks, i € 7*' = {1, 



,nt 



m 



n\AX)^\ql^,,fl^,,{x), X^{x} (1) 

.0, \X\>2 

where q\^, is the probability of existence for the track, 
and flu'i^) is the existence-conditioned PDF 

The full target distribution can be reconstituted from these 
components as |lj p386] 



ft\AX)= J2ft\t'iY)f!iAx-Y) 



(2) 
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where f".^,{X) is the Poisson distribution of undetected tar- 
getsgy[jp366] 
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and /ff ,, {X) is the multi-Bernoulli distribution of targets under 



it\t 



track: 

fP^A{xu...,Xn})^J2 n /*|t'(^-w) (4) 

where the sum is over all functions a : Tf -> {0, . . . , n} 
such that {1, . . . ,n} C a{Tt') and if a{i) > 0, i ^ j then 
a{i) 7^ a.{j), and 



X, 



a{i) 



A ){Xa{€}}, a(«) > 

0, a{i) = 



The PPP intensity of undetected targets is predicted in the 
same manner as a PHD filter ||2l4l 



(5) 
and updated in the same way as a PHD if there are no 
measurements present 



A^|,(x) = [l-P,'>(x)]A?|,_i(x) 



(6) 



We only describe the aspects of the multi-Bernoulli com- 
ponent of the distribution f^.^,{X) that relate to the inter- 
action with the undetected target density; a full derivation 
can be found in |4|. Under the marginal track filter (MTF), 
a new multi-Bernoulli component i is introduced for each 
measurement z in each scan. Using the hypothesis 5 that 
the measurement z is associated with a pre-existing track, 
the hypothesis-conditioned distribution is flu{^) = 1' '•^•' 
probability of existence is zero. Under the hypothesis a that 



'We leave off conditioning on the measurement set history Z* = 
{Zi, . . . , Z^i) tliroughout, as this is implicit in the second subscript ftw- 
^i.e., tai'gets that have never previously been detected 



the measurement is not associated with any previous track, the 
hypothesis-conditioned distribution is Bernoulli, with 



^t\t 






ft\t (^) ^ 
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(7) 
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where 






i[Mz\-)Pi] ^ I ftiz\x)Piix)X'',^,_,{x)dx 



is the Poisson intensity of measurements arriving from previ- 
ously undetected targets. The MTF operates by calculating the 
marginal probability of each measurement-track association 
event (considering aU joint association hypotheses), and rep- 
resenting the posterior distribution via Bernoulli components, 
where each is averaged over the marginal association events 
for that track. In the case of interest, this will simply result in 



Ql\t = p' ia)ql\t ft\t(x) = ft\t (^) 



(9) 



where p*(a) is the marginal probability that the measurement 
z is not associated with any pre-existing track (p'(a)+p*(a) — 
1). As described in [4|, this marginal probability incorporates 
the weight factor w^j^ = Xf{z) + X''^^^_^[ft{z\-)P^]; this 
factor replaces the more common false alarm intensity Xf{z) 
in the standard JPDA association probability calculation. As 
described in [41, we utilise an efficient approximation of 
the joint association probabilities, calculated via loopy belief 
propagation (LBP) (i^. 

In summary, the interaction between the undetected target 
component and the multi-Bernoulli, previously-detected target 
component is 

• The probability of existence of a new Bernoulli compo- 
nent started on a measurement is calculated as the ratio 
between the new target measurement intensity and the 
false alarm measurement intensity plus the new target 
measurement intensity (weighted by the association prob- 
ability). Thus the initial probability of existence will be 
low if the new target intensity is much lower than the 
false alarm intensity, and high in the opposite case. 

• The kinematic distribution of a new Bernoulli component 
started on a measurement is proportional to the unde- 
tected target intensity multiplied by the state-dependent 
detection probability and the measurement likelihood. 
Thus a method is provided for incorporating prior in- 
formation in the track initialisation. 

• In calculation of the marginal association probability, the 
false alarm intensity in the standard JPDA expression is 
replaced by the sum of the false alarm intensity and 
the new target measurement intensity. Accordingly, a 
measurement in an area of high new target intensity is less 
likely to be associated with a previously existing track. 
This provides a principled approach for initiation of new 
targets in the vicinity of existing tracks. 



III. Modelling undetected targets 

In this section, we explore two simple cases in which there 
is an advantage associated with modelling a distribution of 
undetected targets. To commence, we consider the simplest 
case, with uniform birth, death and probability of detection. 
We subsequently consider an adaptive, non-homogeneous case, 
which is representative of a UAV sensor. 

A. Example: Improved initialisation in uniform/stationary 
case 

The first problem that we consider is the most common case 
studied in academia, in which the probability of detection, 
birth probability and death probability are all uniform over 
the region of interest and constant with time. Even under these 
assumptions, we show that it is advantageous to represent a 
distribution of undetected targets, which can be conveniently 
parameterised via a single scalar value. 

The advantage comes from the initialisation of the tracker. 
The most common initialisation does not incorporate any 
tracks, and effectively assumes that no targets are present when 
the system commences operation, and that targets gradually 
arrive from that time in accordance with the birth intensity. 
However, the fact that the tracker has not been operating 
generally does not imply that targets are not present. 

The Poisson distribution of undetected targets provides a 
method for incorporating prior information on the expected 
number of targets present, which is taken into account in calcu- 
lating the existence probability of new Bernoulli components 
as new measurements are received. Perhaps the most natural 
initialisation of the distribution of undetected targets is the 
steady state distribution when the sensor is not active. Under 
the uniform and constant in time assumptions, this is found 
by repeatedly applying the prediction equationr] 






A'' + 



P'K\t^ 



Clearly, this will eventually reach steady state with A"^ = j-rpr 
Thus ApiQ — A"s is the natural initialisation of the tracker (with 
To = 0, i-e., no multi-Bernoulli components). 

We illustrate the advantage that this provides through a 
simple scenario in which targets arrive uniformly in the 
region [—100, 100]^ and velocity region [—1, 1]^, such that 
the hypervolume of the state space is \X\ — 200^ x 2^. The 
state is x = [px^ Vx, Py, Vy]'^ . We set the birth intensity 
to A*" = COS/jA"! (i.e., on average a new target arrives every 
20 time steps) and the survival probability to P^ = 0.999 
(i.e., on average targets survive 1000 time steps; targets and 
hypotheses are removed when they depart the position region 
[-100, 100]2). Thus XI = 50/\X\ (i.e., the expected number 
of targets present is 50). We study a problem involving a 
low probability of detection, P'^ — 0.3, such that several 
measurements are required in order to accrue confidence in the 

^We also require the transition kernel to satisfy (approximately, at least) 
fx ft\t-l{^W)dx' = 1. 



(a) Evolution of expected number of undetected targets 
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(b) Impact of undetected target intensity on tracking performance 
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Fig. 1 . Performance benefit achieved by maintain undetected target intensity. 
Upper figure (a) shows evolution of total expected number of undetected tar- 
gets over time, commencing from the steady state value (50) when the sensor 
is not operating to the steady state value 0.1663 when the sensor is operating. 
Lower figure shows the tracking performance (measured by MOSPA) achieved 
by dynamically estimating the expected number of undetected targets, versus 
assuming different fixed intensities. 



probability of existence of a track. Measurements of position 
are received, i.e., p{z\x) — J\f{z; Hx, R}, where R — I and 

10 
10 



H 



We consider the measurement space to be Z = [—100, 100]'^. 
The false alarm intensity is A*^^ = 10/|-Z| on Z 
(i.e., the expected number of false alarms per scan is 
10)|j For the calculation of the updated kinematic dis- 
tribution of a track upon initialisation, we approximate 
K\t-ii^) ~ iK\t^i ■ \X\W{x;0,P} where P = 
diag[100V3, 1/12, 100V3, 1/12]. This approximation 
matches the mean and covariance of the uniform distribution 
on X. the target kinematics follow the widely-used model with 
velocity as a random walk, with diffusion q = 0.01. The multi- 
Bernoulli tracker is described in |4|; it outputs all tracks with 
a probability of existence of at least 0.8, and for which the 
trace of the covariance is less than 10. 

The simulation runs for 100 time steps. The expected 
number of undetected targets is shown in Fig. [TJa) as a 
function of time. As expected, the initial number of undetected 
targets decays from a large value initially, to steady state 
operation by approximately t = 30. In Fig. flTb), we show the 
average performance as measured by the mean optimal sub- 
pattern assignment (MOSPA) metric 1 10] (measured against 
the full state vector), with p = 2 and c = 10. We compare 
the proposed dynamic method for estimation of the undetected 
target intensity to alternatives which hold the intensity of unde- 
tected targets constant at X'^u_i ■ \X\ = 5, 1, and 0.1663 (the 

'^The measurement space is chosen to be larger than the birth region in 
order to avoid targets departing the portion of the space where false alarms 
exist. The tracker assumes that both the undetected target intensity and false 
alarm intensity are constant throughout the space. 



Steady state value with the measurement process operating). 
The figure shows that the method incorporating the dynamic 
estimate of undetected target intensity outperforms the steady 
state alternatives throughout the simulation. As expected the 
filter assuming the steady state intensity performs well at 
the end of the simulation (when the steady state condition 
prevails), whereas the filters assuming higher intensity perform 
better (than other steady state alternatives) during the initial 
transient. Dynamically modelling the expected number of 
undetected targets can speed the process of initially acquiring 
tracks in a scene without trading off steady state performance. 

B. Example: Adaptive initiation for moving sensor 

In the more general (and rarely studied) case in which the 
distribution of undetected targets is non-homogeneous, an 
appropriate representation must be chosen. We propose dis- 
cretising the state space in order to represent this distribution 
via a grid filter. Grid-based representations (e.g., [llj ) have 
fallen out of vogue for most estimation applications as they are 
generally inefficient. In particular, in any problem where the 
probability distributions are peaked, grid representations result 
in large computing resources being spent on the vast majority 
of cells with near-zero probability density. In comparison, the 
preferred sample-based methods allow computing resources to 
be focussed on the region with significant probability mass. 
However, in the present application, the Poisson intensity 
function of undetected targets is, by nature, both diffuse and 
smooth. Thus a grid-based representation is efficient (since 
most grid points will have comparable intensity), and, farther- 
more, a coarse discretisation suffices, such that the grid-based 
method is tractable. Conversely, sample-based methods using 
any reasonable number of samples would be likely to exhibit 
poor sample coverage in the vicinity of many cells. 

To demonstrate the more general case, we consider a 
problem in which the sensor is moving, and the probability 
of detection is zero outside of a cone 45° either side of the 
sensor's heading. The sensor trajectory is shown in Fig. l2] in 
blue (it commences at (—20, —20), heading directly down the 
map). All other parameters remain identical to Section |III-A| 
except that there are two grid cells (shown as magenta stars 
in Fig. l2| in which the birth rate is equivalent to 20% 
of the birth rate in the entire region in the previous case 
(e.g., representing airports). False alarms in regions where the 
detection probability is zero are naturally discarded by the 
update equations. 

The discretisation used to represent the intensity of un- 
detected targets centres cells on {—100, —96, ... ,96, 100} 
in position dimensions, and {-1,-0.6,-0.2,0.2,0.6,1} in 
velocity dimensions, for a total number around 94, 000 cells. 
The transition kernel is obtained by offline Monte Carlo 
simulation, drawing large number of samples uniformly within 
the cells, simulating the continuous dynamics process, and 
observing in which cell they arrive. The survival probability 
and birth density are modified at the region boundary to 
achieve a uniform steady state intensity (in the absence of 
the "airport" entry points in the birth intensity). Fig. l2] shows 
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Fig. 2. Undetected target intensity at time t = 110 (darker shades denote 
higher undetected target intensity). Sensor trajectory is show in blue, and field 
of view at t = 110 is shown in green (P"* = 0.3 within the field of view, 
and zero outside). 
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Fig. 3. Performance benefit obtained by maintaining an adaptive, non- 
homogeneous estimate of undetected target intensity for case of a moving, 
manoeuvring sensor. Blue line shows performance with dynamic estimate, 
versus alternatives assuming different fixed, uniform intensities. 



a snapshot of the scenario, at time t — 110 (the scenario 
duration is 140 time steps). The background colour of the 
figure shows the undetected target intensity (summing over ve- 
locity dimensions); lighter shades denote lower intensity, while 
darker shades denote higher intensity. The region presently 
being observed has the lowest intensity while the area that 
has not been observed for the longest duration has the highest 
intensity (as expected). Additionally, intensity is increased near 
the region boundary, modelling entry of new objects into the 
region. 

The tracking performance for the scenario is shown in 
Fig. [3] The diagram again shows the MOSPA measure of 
tracking performance, this time counting only those targets 
that are within the sensor coverage region. The diagram shows 
that, each time the sensor manoeuvres, a performance transient 
occurs which is similar to resetting the tracker The proposed 



method which maintains a dynamic, non-homogeneous esti- 
mate of undetected target intensity is exhibits an significant 
improvement in its ability to adapt to these transients. Again, 
the consequence of this is faster track initiation during this 
transient, without compromising the rate of false track occur- 
rence at other times. 

IV. Recycling: Hybrid PHD-MeMBer 

A. Derivation of method 

Recall again (from Section II-B i the form that the MTF 
maintains from one time interval to the next. In probability 
generating functional (PGFl) LL p371] form, the overall tar- 
get distribution (incorporating both undetected and previous 
detected targets) is Q 



(10) 



where 



K\M 



^(2;)Aj|t(a;)da:: 



The PGFl form makes clear that this is a union of independent 
components: a Poisson point process component, and a series 
of independent multi-target Bernoulli componentsP] Since a 
new track must be created for every new measurement re- 
ceived, tracks will inevitably need to be deleted from the 
system. The system designer must trade off system perfor- 
mance against computational complexity in constructing this 
mechanism. The major source of low probability of existence 
tracks is initiation. If the probability of detection is low and the 
false alarm rate is high, then many tracks will be created on 
false alarms, and it will take a large number of measurement 
scans to reduce the probability of existence of the tracks to 
the point where they can be excluded with confidence. 

One possibility which arises from the form of ( fTO] ) is to 
approximate individual tracks as being Poisson, rather than 
deleting them. Consider a single Bernoulli component ]\^{X) 
of the form (m), parameterised by q\^ and f\^{x). We may 
choose to approximate this component as being Poisson: 

/,»l,(X)=exp(-AJ|jn^t|*(^) 



x^X 



where A* 



y^U — /AJ|j(2;)dx, and we define Jluix) = 
Aj|j(x)/AJ|j for later use. The distortion caused by this approx- 
imation may be measured by the multi-target KL divergence: 
|l,p513] 



D{nu{x)\\n\t{^)) 



(i-<?jit)iog 






it\t 



exp(-AJ|J 



The multi-target Bernoulli components are not naturally independent but 
they are approximated as such |4|, as in the MeMBer il|3i. 



KL divergence between Bernoulli and Poisson 
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Fig. 4. Multi-target KL divergence between multi-target Bernoulli dis- 
tribution and best-fit Poisson distribution as a function of target existence 
probability g*,^. 



ql\tft\ti^)^o&- 



il\tfi\ti^) 



-dx 



= (l-gJ|,)log(l-gj|,) + (l-gJ|JAJ|,+ 

+9J|Jlogg,V + Aj|* - log Aj|,] + <zj|,i?(/;|,(x)||/,^|,(x)) 

As shown in fl] p579], it is optimal to set 

4(x) = i?|,(:r), AJ|, = gj|„ Xl^,{x) ^ ql^Jl^.ix) (11) 

The value of the KL divergence at this optimal choice is: 

D{fl^,{X)\\fl^,{X)) = gj|, + (1 - gj|,) log(l - gj|J 

The value of this KL divergence is shown in Fig. |4] as a 
function of the target existence probability q\^. The figure 
demonstrates that the distortion caused by the approximation 
(as measured by the KL divergence) is very small for exis- 
tence probabilities less than 0.1 (when q. 



t\t 



0.1 we have 



D{fl\ti^)\\ft\ti^)) ~ 0.005, and when ql^^ = 0.2 we have 
^(/*itWll4W)«0.02). 

Theorem [l] shows that the KL divergence between the 
overall multi-target distribution comprised of independent 



components (e.g., ( 10 1) and a modified multi-target distribution 
in which approximations have been made to a number of these 
components is bounded by the sum of the KL divergences 
between the components and their respective approximations r] 

Theorem 1. Let 



fix) = ^ g{W)h{X - W) 
w<zx 

fix) = Y. ~9(W)hiX - W) 
wcx 

Then DifWf) < Dig\\~g) + D{h\\h). 

The proof of this theorem is in the appendix. Thus the over- 
all distortion we cause to the complete multi-target distribution 
by approximating a number of components is bounded by the 
sum of the individual component distortions, which depend 
only on the component existence probabilities. Accordingly, 

*The theorem is stated and proven for the two component case; the general 
case follows simply by induction. 



we may approximate the components with lowest existence 
probabilities such that the sum of the distortions is less than 
an overall distortion budget. 

When any number of Bernoulli tracks are approximated 
as being Poisson, the resulting multi-target distribution is 
equivalent to one in which those tracks are dropped, and 
their intensity is added onto the undetected target intensity. 
To confirm this, denote the subset of tracks that we retain as 
Tt; the approximated distribution is then 

G,|,[/i](xexp{A?|J/.]}. n GJ|J/i] • n Gj|j/i] 

'l£Tt\ft i€ft 

(xexp{A?|J/i]}. n exp{AJ|J/j]} . n G?|,[/i] 
leTtXft ieft 



exp< A^iJ/i] 



E 

ieTt\ft 



H\t 



m-X{G\\M 



ieTt 



expW|J/i]}-nGJ|,W 



where 



A"|*N = A^|,W 
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ieTt\ft 



K\t[h] 



or, in intensity form. 



E 

ieTt\ft 



Kiti^) 



We refer to this concept as recycling, since the tracks that we 
delete are re-used in generation of new tracks in subsequent 
measurement scans. Not surprisingly, it can be shown that if 
the prior distribution is purely Poisson (i.e., there are no pre- 
existing target tracks) and we choose to recycle all posterior 
tracks, then the posterior Poisson distribution is equivalent 
to that obtained using the PHD. This can be easily seen by 
observing that, if there are no prior tracks, p*(a) = 1, hence 
by (|9) and (fTT), 



ftiz\x)P,^ix)\-^^_^ix) 
A'/(z) + A^|,_J/t(z|.)Pf 



K\ti^) 



Combining this with |6|) and setting 7t = (i.e., recycling 
all tracks), we obtain an update for a measurement set Zt — 



Ul-- 



'} of 



A^|,(x) = [l-P,''(x)]A^I,„,(x)- 



2^ \fa 



ft{zi\x)PHx)\'iu,ix) 



^^^\f{zi) + \''^^^^^[f,izi\-)P,'] 

which is the PHD update f2l. 

By recycling a subset of tracks, we permit the large mass 
of tracks with low probability of existence to be represented 
efficiently by the Poisson distribution, while maintaining ex- 
plicit Bernoulli tracks on the subset with non-negligible prob- 
ability of detection. Representing low probability of existence 
tracks via the Poisson distribution reduces the computational 
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Fig. 5. Performance of recycling-based method, compared to alternatives. 
Top figure shows tracking performance measured by MOSPA, while bottom 
figure shows average number of tracks maintained. 



burden due to data association. Furthermore, if the Poisson 
distribution is represented as a discrete grid, then there is 
no computational cost associated with representing additional 
tracks. 

In practice, this approximation allows the system to grad- 
ually accrue confidence in the presence of a target before 
choosing to maintain an explicit track. As shown in (|7]i, the 
existence probability of a new, isolated track is the ratio 
between the undetected target intensity (in the vicinity of 
the measurement) and the undetected target intensity plus the 
false alarm intensity, so the intensity added by recycling will 
cause the existence probability of a new track due to a later 
measurement in the same vicinity to be increased, reducing 
the likelihood that the track will again be recycled. 

B. Example: Reduction of tentative tracks 

We demonstrate the utility of recycling using the same scenario 



as Section III-A and the same birth intensity and discretisation 
as Section III-B Following the result in Fig. HI we project 
tracks with a probability of existence q^.^ < 10~^ onto the 
Poisson component (i.e., we recycle them). We compare the 
performance of the method utilising recycling to alternatives 
which delete tracks with (/L^ < 77, where 77 — 10^^ or 77 = 
10^^ (the latter value was used in the simulations in Sections 
[HTaI and [II^B]). 

Fig. I5] demonstrates the utility of the recycling. The top 
diagram shows the tracking performance, measured by the 
MOSPA metric. The method utilising recycling achieves simi- 
lar (and even slightly improved) performance to the pure multi- 
Bernoulli method with rj = 10^"^. The performance of the 
multi-Bernoulli method with 77 = 10^^ reduces through the 
simulation as the undetected target intensity lowers, and the 



tracker becomes unable to initiate tracks on new targets (as 
they are deleted before being able to be confirmed). 

The practical advantage of recycling is in the number 
of multi-Bernoulli components that need to be maintained. 
Fig.|5|b) shows the average number of components maintained 
by each method versus time. The method utilising recycling 
maintains slightly more tracks than the method with 77 = 10^^, 
or around a quarter of the tracks of the method with 77 = 10^'^. 
Thus recycling permits similar performance to the alternative 
with a very low track deletion threshold, using only a fraction 
of the number of tracks (and, in turn, a reduction in compu- 
tation for gating, data association, etc). 

The result is somewhat surprising since the discretisation 
used for the Poisson component was quite coarse, e.g., the 
measurement variance is 1, yet the discretisation resolution 
is 4. Thus the method appears to be quite robust to the 
discretisation used in a grid representation. 

Finally, we note that the grid representation is not necessary 
for the algorithm; indeed, a sample-based representation is 
also possible (as in [Tl), as is a hybrid approximation (e.g., 
retaining Gaussian mixture components until their variance 
grows sufficiently to be adequately represented by a grid). 

V. Related work 



The landmark paper |12| incorporated a Poisson distribution 
of undetected targets. Its primary role was to provide a 
distribution of newly detected targets for multiple sensors with 
overlapping fields of view, and as such, it was approximated 
as being uniform at the beginning of each scan|j The RFS 
framework [IJ has produced an improved set of analytical 
tools for Bayesian modelling in tracking problems, and the 
birth, transition and death models we employ are standard in 
the PHD. Despite this resurgence, the concept of maintain- 
ing a distribution of undetected targets appears not to have 
reentered the mainstream. The Gaussian mixture and sample- 
based representations commonly used to represent the PHD 
generally struggle to adequately model the diffuse distribution 
of undetected targets. 

To the author's knowledge, the concept of mixing multi- 
Bernoulli and PHD representations was first proposed simul- 
taneously in the pre-print |6| and in |7l. The latter work 
was somewhat heuristic, in that it was not supported by the 
RFS-based derivation in | |4|6| . This, for example, resulted in 
retaining the standard hypothesis weight for a target not being 
associated to any other measurement, e.g., w^f^ = Xf{z), as 
opposed to the value in the previous section, which results 
from the RFS derivation. The analytical tools provided by the 
RFS framework also permit an analysis of the error cause by 
a Poisson approximation to a multi-Bernoulli component (i.e.. 
Fig. Hll, which permits a system designer to set the threshold 
for when to use the approximation. 

A significant contribution in fj] was the development of 
a sophisticated sample-based method for representing the 
PHD intensity. A surprising outcome of the experiment in 

^Presumably, this is due to the computational limitations of the day. 



Section IV-B is that a grid-based approximation of the intensity 
with a resolution several times lower than the sensor accuracy 
still appears to be adequate to achieve equivalent performance 
to a system using a very low track deletion threshold. 

VI. Conclusion 

This paper has shown the practical advantage of hybrid 
tracking systems that utilise both a Poisson component and 
a multi-Bernoulli component. Our sister paper |4| showed 
that this is a natural result of the derivation of the Bayes 
RFS filter, and that the Poisson component represents the 
intensity of targets that remain undetected. Here we have 
demonstrated the practical advantage of maintaining such a 
quantity, improving track initiation performance in dynamic 
environments. Subsequently, we derived the method of re- 
cycling, which additionally utilises the Poisson component 
to represent Bernoulli components with low probability of 
existence, and demonstrated its practical utility in tracking 
problems. 

Appendix 
sub-additivity of kl divergence 

This section presents the proof of Theorem [T] which in turn 
requires Lemma [T] 

Proof of Theorem 17} 
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where (a) is a consequence of the log-sum inequality, 1 13 
p29], and (6) is the result of Lemma [T] I 
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Proof: For simplicity, when n = 
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where (a) results from the observation that terms in the sum 
over / with the same number of elements in / will have the 
same value, and (6) results from letting I = [n ~ m) and 
reordering terms in the two-dimensional summation. ■ 
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